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Abstract 

The dark energy models with variable equation of 
state parameter co are investigated by using law 
of variation of Hubble 's parameter that yields the 
constant value of deceleration parameter. Here the 
equation of state parameter uj is found to be time 
dependent and its existing range for this model is 
consistent with the recent observations of SN la 
data, SN la data (with CMBR anisotropy) and 
galaxy clustering statistics. The physical significance 
of the dark energy models have also been discussed. 
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1 Introduction 

Observations on large scale structure (LSS) [T] and 
cosmic microwave background radiation [2] — [5] indi- 
cate that the Universe is highly homogeneous and 
isotropic on large scales. Under standard assump- 
tions on the matter content, in general relativistic 
cosmological models, isotropy is a special feature, re- 
quiring a high degree of fine tuning in order to repro- 
duce the observed Universe. Recent years have wit- 
nessed the emergence of the idea of an accelerating 
Universe due to some observational results [7] — [10]. 
This signifies a remarkable shift in cosmological re- 
search from expanding Universe to accelerated ex- 
panding Universe. Now, the problem lies in detect- 
ing an exotic type of unknown repulsive force, termed 
as dark energy (DE) which is responsible for the ac- 
celerating phase of the Universe. The detection of 
DE would be a new clue to an old puzzle: the grav- 
itational effect of the zero-point energies of particles 
and fields [11]. The total with other energies, that 
are close to homogeneous and nearly independent of 
time, acts as DE. The paramount characteristic of 
the DE is a constant or slightly changing energy den- 
sity as the Universe expands, but we do not know 
the nature of DE very well (see [12] — [2^ for reviews 
on DE) . DE has been conventionally characterized by 
the equation of state (EoS) parameter lu = p/p which 
is not necessarily constant. The simplest DE candi- 
date is the vacuum energy (w = —1), which is argued 
to be equivalent to the cosmological constant (A) 
[21j . The other conventional alternatives, which can 
be described by minimally coupled scalar fields, are 
quintessence (w > —1), phantom energy (io < — 1) 
and quintom (that can across from phantom region 
to quintessence region) as evolved and have time 
dependent EoS parameter. Some other limits ob- 
tained from observational results coming from SN la 
data [22] and SN la data collaborated with CMBR 
anisotropy and galaxy clustering statistics [23j are 
-1.67 < w < -0.62 and -1.33 < w < -0.79 respec- 
tively. However, it is not at all obligatory to use a 
constant value of uj. Due to lack of observational ev- 
idence in making a distinction between constant and 
variable w, usually the equation of state parameter 
is considered as a constant [M] [35] with phase wise 
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value —1, 0, 1/3 and +1 for vacuum fluid, dust fluid, 
radiation and stiff fluid dominated universe, respec- 
tively. But in general, w is a function of time or 
redshift [IS] — [IH]- For instance, quintessence models 
involving scalar fields give rise to time dependent EoS 
parameter w [5^ — [22]. Also some literature is avail- 
able on models with varying fields, such as cosmolog- 
ical model with variable equation of state parameter 
in Kaluza-Klein metric and wormholes [33] [31]. In 
recent years various form of time dependent w have 
been used for variable A models [551 13S]- Recently 
Ray et al [37], Akarsu and Kilinc [38] have studied 
variable EoS parameter for generalized dark energy 
model. 

Bianchi type-V universe is generalization of the 
open universe in FRW cosmology and hence its 
study is important in the study of DE models in a 
universe with non-zero curvature |39| . A number of 
authors such as Collins [ID], Maartens and Nel [H], 
Wrainwright et al [H], Canci et al [33], Pradhan et 
al [33] and Yadav [3S] have studied Bianchi type-V 
model in different physical contexts. Recently Yadav 
and Yadav [46] have studied anisotropic DE models 
with variable EoS parameter. In this paper, we 
have investigated the DE models with variable uj in 
Bianchi type-V Universe. This paper is organized as 
follows: The metric and field equation are presented 
in section 2. In section 3, we deal with the solution 
of field equations and discussion. Finally the result 
are discussed in section 4. 



2 The Metric and Field Equa- 
tions 

We consider LRS Bianchi type V metric in the form 

ds^ = ~dt^ + A^dx^ + B^e^'-" {dy^ + dz^) (1) 

where A and B are the function of t only. 
The simplest generalisation of EoS parameter of per- 
fect fluid may be to determine the EoS parameter 
separately on each spatial axis by preserving the di- 
agonal form of the energy-momentum tensor in a con- 
sistent way with the considered metric. Thus, the 



energy momentum tensor of fluid is taken as 

Tl =diag[T^,TlTlTi] 
Then, one may parametrize it as follows, 
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diag [p, -px, ~Py, -Pz] = diag [1, -uj^, -^^y, -w^] p 
= diag [1, -w, ~ {lu + S) , - {uj + S)] p (3) 



where p is the energy density of fluid, Px, Py and 
Pz are the pressures and uJx, ijJy and lOz are the di- 
rectional EoS parameters along the x, y and z axes 
respectively, uj is the derivation-free EoS parame- 
ter of the fluid. We have parametrized the deviation 
from isotropy by setting ujx = uJ and then introduc- 
ing skewness parameter 5 that are deviation from uj 
along y and z axis respectively. 

The Einstein field equations, in gravitational units 
(SttG = 1 and c = 1), are 



(4) 



where the symbols have their usual meaning. 
The Einstein's field equation (j31) for the Bianchi type 
V space-time (|T]), in case of (|3]), lead to the following 
system of equations 
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= -{uj + 6)p (7) 
(8) 







Here, the sub indices 4 in A, B, and elsewhere denote 
differentiation with respect to t. 
Integrating equation ([8]), we obtain 



A^kB 
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where k is the positive constant of integration. We 
substitute the value of equation ^ in equation ([7]) 
and subtract the result from equation ([6]), we obtain 
that the skewness parameter on z-axis is null i.e. 

5^0 
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Thus system of equations from ([S]) - may be re- 
duce to 

^— (10) 



B 
' B 



B^ 



1 



52 



-up 



(11) 



Now we have three Hnearly independent equations 
([9]) - (fTTj) and four unknown parameters {A, B, u, p). 
Thus one extra condition is needed to solve the sys- 
tem completely. To do that, we have used the law 
of variation of Hubble's parameter that yields a con- 
stant value of deceleration parameter. 
The average scale factor of Bianchi type V metric is 
given by 

R={AB^e^'')^ (12) 

We define, the generalised mean Hubble's parameter 
H as 

H^^iHi+H2+ H3) (13) 

where Hi = H2 = H3 = ^ are the directional 
Hubble's parameter in the direction of x, y and z 
respectively. 

Equation (13), may be reduces to 
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Here the line-element ([T]) is completely characterized 
by Hubble's parameter H. Therefore, let us consider 
that mean Hubble parameter H is related to average 
scale factor R by following relation 



H = kiR- 



(15) 



where ki > and n > 0, are constant. Such type 
of relation have already been considered by Herman 
[U] for solving FRW models. Later on many authors 
(Singh et al |48n49j and references therein) have stud- 
ied flat FRW and Bianchi type models by using the 
special law of Hubble parameter that yields constant 
value of deceleration parameter. 
The deceleration parameter is defined as 
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From equations and (15), we get 

R4 = kiR-'^+^ (17) 
Ra4 = -kl{n - l)i?-2"+i (18) 
Using equations (IT4l) and (ITS)) , equation (fT6)) leads to 



q = n — I (constant) 



(19) 



The sign of q indicates whether the model inflates or 
not. The positive sign of q corresponds to standard 
decelerating model where as the negative sign of q 
indicates inflation. It is remarkable to mention here 
that though the current observations of SN la (Per- 
mutter et al 1999, Riess et al 1998) and CMBR favour 
accelerating models i.e. q < 0. But both do not 
altogether rule out the decelerating ones which are 
also consistent with these observations (Vishwakarma 
2003). 

From equation (jl7p . we obtain the law of average 
scale factor R as 
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{Dt + ci) 



when n ^ 
when n = 



(20) 



where Ci and C2 are the constant of integration. 
From equation ([20]) . for n 7^ 0, it is clear that the 
condition for accelerating expansion of Universe is 
< n < 1. 



3 Solutions of the Field Equa- 
tions and Discussion 

3.1 Case (i): when n ^ 

Equations dH), ^ and ([20]) lead to 



(21) 



B ^l{Dt + ci)"- 
From equations (jH]) and (|21l) . we obtain 

A ^ L{Dt + ci)" 

where 1 is the constant of integration and L — kl. 
Thus the Hubble's parameter (iJ), scalar of expan- 
sion (0), shear scalar (cr^) and spatial volume [V) are 
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given by 



H = 



ki 



(Dt + ci) 
3fci 



3fc2 - 
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(23) 
(24) 

(25) 
(26) 



Using equations ([TOl) . (f20| and ([2T|) . the energy den- 
sity of the fluid is obtained as 



(27) 



Using equations (HH), ^ and (gT]), the equa- 

tion of state parameter cj is obtained as 



(L't + ci)" 



2n 



(IJt + Cl) 



(28) 



From equation ((27|) . we note that p(t) is the 
decreasing function of time. This behaviour is 
clearly shown in Fig. 2, as a representative case 
with appropriate choice of constants of integration 
and other physical parameters using reasonably well 
known situations. From equation (j28p . it is observed 
that the equation of state parameter w is time 
dependent, it can be function of redshift z or scale 
factor R as well. The redshift dependence of oj can 
be linear like uj{z) = ujo + w^z with uj^ = {^) = 
[501 EI] or nonlinear as w(z) = ujq + feflSS]. 
The SN la data suggests that -1.67 < w < -0.62 
while the limit imposed on by a combination 
of SN la data (with CMB anisotropy) and galaxy 
clustering statistics is —1.33 < w < —0.79 So, 
if the present work is compared with experimental 
results mentioned above then, one can conclude 
that the limit of uj provided by equation psp may 
accommodated with the acceptable range of EoS 
parameter. The value cj = — 1 is the case of vacuum 
fluid dominated Universe. But here, we are dealing 
with the solution for rt 7^ 0, therefore in this case, 
vacuum fluid dominated universe is meaningless. 




Figure 1: The plot of spatial volume V for n = .8 



Also we see that during evolution of Universe, at an 

— Cl ) , the 



instant of time t ^ j) 



W2 



w vanishes. Thus at this particular time, our model 
represents dusty Universe. 

The critical density (pc) and density parameter (fl) 
are given by 

3fc2 



D^L^ -n^Dt + ciY 



(29) 



(30) 



The variation of equation of state parameter (w) with 
cosmic time (t) is depicted in figures 3 and 4 as a 
representative case with appropriate choice of con- 
stants. Fig. 3 and 4, clearly show that w is evolving 
with negative value and the existing range of w is 
in nice agreement with SN la data [32]. Thus our 
model is a realistic model. The plots of uj for n=0.5 
and n=0.85 indicate that u merge well with SN la 
and CMBR observations [Fig. 3 and Fig. 4]. The 
variation of density parameter (51) with time in ac- 
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Figure 2: The plot of energy density (p) vs. time (t) 



Figure 4: The plot of EoS parameter (cj) vs. time (t) 
for n = 0.85 
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Figure 6: The plot of EoS parameter (uj) vs. time (t) 
for n = 



V = Lo/^e(3fci*+2-) (36) 

Using equations pO)) . (|3T|) and ([32]) . the energy den- 
sity of the fluid is obtained as 



P = 



,-2kit\ 



7^2 
^0 



(37) 



Using equations ([TT|) . (|3T|) . ([32]) and ([37]). the equa- 
tion of state parameter uj is obtained as 



-2A;ii 



3(fc2L2_e-2'=i*) 



(38) 



The variation of the EoS parameter uj with cosmic 
time (t) is shown in Fig. 6. The value of u is found 
to be negative which is supported by SN la data and 
galaxy clustering statistics [22l |23] . 
The critical density (pc) and density parameter (fl) 
are given by 

Pc = 3kl (39) 



n = 



klL^ — e 



-2feit 



(40) 



celerating mode of Universe is clearly shown in Fig. 
5. 

3.2 Case(ii): when n = 

Equations ©, (15) and ^ lead to 

B = loe'''* (31) 

From equations (|9]) and (|3T|) . we obtain 

A ^ ioe''^* (32) 

where Iq is constant of integration and Lq = klo- 
Thus the Hubble's parameter {H), scalar of expan- 
sion (9), shear scalar (ct^), and spatial volume {V) 
are given by 

H^ki (33) 
6 = 3/ci (34) 
cr^ = (35) 



4 Concluding Remarks 

In this paper, we have studied dark energy model 
with variable EoS parameter uj. The model is derived 
by using law of variation of Hubble's parameter that 
yields a constant value of deceleration parameter. 
We have studied two cases, (3.1) and (3.2) for n ^ 
and n = respectively. In both the cases, uj is found 
to be time varying and negative which is consistent 
with recent observations (22] [23]. In case (i), we 
have shown that, lu is evolving from uj < —1 and 
finally end up with > — 1, representing two phases 
of universe i.e. a; < — 1 (Phantom fluid dominated 
universe) and uj > —I (quintessence), [Fig. 3 and Fig. 
4]. Where as in case (ii), for n = 0, from equation 
([55]) . we have obtained that, at cosmic time t = 0, 
UJ < —1 and when t — ^ oo, uj — >■ —1. Therefore, 
one can conclude that at early stage, universe was 
dominated by phantom fluid and at late time it will 
be vacuum fluid dominated universe [Fig. 6]. Also 
it is seen that that in both cases, ^ — constant, 
therefore the proposed models do not approach to 
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isotropy at any time. 

Here the age of universe is given by 




which is also differ from the present estimate i.e. 
To = Hq^ « 14Gyr. But if we take, /ci = 1 then our 
model {n ^ 0) is in good agreement with present age 
of universe. 

The main fetures of the models are as follows 

• Though there are many candidates such as 
cosmological constant, vacuum energy, space-time 
curvature, cosmological nuclear energy, etc as re- 
ported in the vast literature for DE, the proposed 
models in this paper at least present a new candidate 
(EoS parameter) as possible suspect of the dark 
energy. 

• The dark energy models are based on exact so- 
lution of Einstien's field equations for the Bianchi 
type V space-time filled with perfect fluid. To our 
knowledge, the literature hardly witnessed this sort 
of exact solution for Bianchi type V space-time. So 
the derived DE models add one more feather to the 
literature. The DE models presents the dynamics 
of EoS parameter uj provided by equations p8| and 
((38)) may accommodated with the acceptable range 
-1.67 <Lu < -0.62 of SN la data (Knop et al 2003). 
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